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We present a controlled method for computing the exchange coupling in correlated one- 
dimensional electron systems based on the relation between the exchange constant and the pair- 
correlation function of spinless electrons. This relation is valid in several independent asymptotic 
regimes, including low electron density case, under the general condition of a strong spin-charge 
separation. Explicit formulas for the exchange constant are obtained for thin quantum rings and 
wires with realistic Coulomb interactions by calculating the pair-correlation function via a many- 
body instanton approach. A remarkably smooth interpolation between high and low electron density 
results is shown to be possible. These results are applicable to the case of one-dimensional wires of 
intermediate width as well. Our method can be easily generalized to other interaction laws, such as 
the inverse distance squared one of the Calogero-Sutherland-Moser model. We demonstrate excel- 
lent agreement with the known exact results for the latter model and show that they are relevant 
for a realistic experimental setup in which the bare Coulomb interaction is screened by an edge of 
a two-dimensional electron gas. 
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I. INTRODUCTION 



Much interest has been devoted to the spin degree 
of freedom in novel one-dimensional (ID) conductors, 
both of linear shape (carbon nanotubes,i semiconduc- 
tor nanowireSfS conducting molecules^) and recently, of 
a circular one (quantum ring o"^i^i^ ). Parameters of these 
systems, e.g., average distance between the electrons a, 
their total number N , their effective mass m, dielectric 
constant e, effective Bohr radius as = h^t/me^, etc., can 
vary over a broad range or can be tuned experimentally. 
This creates an intriguing opportunity of reaching the 
Wigner-crystal (WC) limitf24& = a/2aB ^ 1, where 
electrons arrange themselves into a nearly regular lat- 
tice. According to numerical simulations,* the ID WC is 
well formed already at > 4. The corresponding elec- 
tron densities are easily achievable but the presence of 
disorder has so far hindered experimental investigations 
of the ID WC regimeiiS. Remarkably, the unwanted dis- 
order can apparently be substantially suppressed in the 
case of carbon nanotubes suspended above a substrate. 
In such devices very large with no immediately obvious 
disorder effects have been recently demonstrated.^*^ Due 
to a finite length of the nanotubes used, they contained 
only a few electrons (< 30) in the > 4 regime. (Such 
systems are commonly referred to as Wigner molecules.) 
These encouraging progress on the experimental side and 
a recent revival of interest to the ID WC on the theoret- 
ical sidsii have motivated us to undertake in the present 
paper a careful investigation of the fundamental energy 
scales that govern the spin dynamics in ID Wigner crys- 
tals and molecules. 

In a strict academic sense, the term "crystal" should 
not be used in ID because the phonon-like vibrations of 
the putative lattice are infrared divergent. It is never- 
theless true that at rg ^ 1 the pair correlation function 




FIG. 1: (Color online) A sketch of the instanton trajectories 
for the six-electron Wigner molecule on a ring. 
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is sharply peaked at integer multiples of a. Here 
electron coordinates. Thus, for the study of electron cor- 
relations in Wigner molecules or in a group of few nearby 
electrons in an infinite wire the WC concept is fully ad- 
equate. Once adopted, this concept implies^^ that the 
low-energy dynamics of electron spins is dominated by 
the nearest-neighbor Heisenberg interaction of strength 
J proportional to an exponentially small probability of 
quantum tunneling under the Coulomb barriers that sep- 
arate the adjacent electrons. Exchanges of more distant 
neighbors are penalized with a much stronger tunneling 
suppression and can therefore be neglected. 

As a result of the exponential smallness of the exchange 
coupling J, the energy scales for orbital and spin dynam- 
ics are drastically different. This strong spin-charge sep- 
aration is expected to cause anomalies in many essential 
electron properties, e.g., balhstic conductanceii of quan- 
tum wires and persistent current of quantum rings. To 
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systematically explore such phenomena a reliable esti- 
mate of the exchange coupling J is needed. This has 
been an important open problem however. The small- 
ness of J makes it difficult to compute even by rather 
advanced computer simulations.^^ Attempts to derive J 
analytically (for the nontrivial case > 2) have so far 
been based on the approximation that neglects all de- 
grees of freedom in the problem except the distance be- 
tween the two interchanging electrons>ii*i^ We call this 
a Frozen Lattice Approximation (FLA). The accuracy of 
the FLA is unclear because it is not justified by any small 
parameter. When a given pair does its exchange, it sets 
all other electrons in motion, too (Fig. Simple di- 
mensional considerations (confirmed by the calculations 
below) show that the maximum tunneling displacements 
of the electrons nearby the first pair reach a fraction of 
a, i.e., they are at most numerically, but not paramet- 
rically, small. In view of the aforementioned and other 
recent experimentali^*i& and theoretical^^'^*'^^ work on 
spin-related effects in ID conductors a controlled calcu- 
lation of J seems a timely goal. It is accomplished in this 
paper, where we treat the spin exchange in a Wigner 
molecule (or a WC) as a true many-body process and 
compute J to the leading order in the large parameter 7% . 
A brief account of this work has been reported in Ref . 12(1. 
The same problem has been independently studied by 
Klironomos et al.^ and wherever they can be compared 
our results agree. 

The paper is organized as follows. The definition of 
the model and our main results are given in Sec. ^ The 
step-by-step derivation for the simplest nontrivial case 
of a three-electron molecule is presented in Sec. lIIII The 
N > 3 case is discussed in Secs. lIVI andlVl An alternative 
route to the derivation of J is outlined in Sec. IVII The 
comparison with the known analytical results for related 
models is done in Sec. IVIII The issues relevant to the 
current experiments are addressed at the end of Sec. IVIII 
and in Sec lVIIll The concluding remarks are in Sec. lIXI 



II. MODEL AND MAIN RESULTS 

We assume that electrons are tightly confined in the 
transverse dimensions on a characteristic lengthscale 
("radius" of the wire) R <^ as- In this case the energy 
separation h^/niR^ of the ID subbands greatly exceeds 
the characteristic Coulomb interaction energy e^/eR. As- 
suming that only the lowest subband is occupied, the 
Hamiltonian can be projected onto the Hilbert space of 
this subband. This turns the problem into a strictly ID 
one, with the Coulomb law replaced by an effective in- 
teraction U{r) that tends to a finite value at distances 
r <C i?. The precise form of U depends on the details of 
the confinement. For simplicity, we adopt a model^SiSi 

e2 1 

U{r)^ — -, (2) 

e r + R 



iV 3 4 6 8 00 00-FLA 

J7 2.8009 2.7988(2) 2.7979(2) 2.7978(2) 2.7978(2) 2.8168 
K 3.0448 3.18(6) 3.26(6) 3.32(7) 3.36(7) 2.20 



TABLE L Results for Wigner molecules on a ring (finite A'') 
and for wires {N = 00). 

which nevertheless correctly captures both short- and 
long-range behavior of the interaction for any realis- 
tic confinement scheme and is similar to other expres- 
sions used in the literature^S^iS^ see more in Sec. IVIIII 
Most of our calculations below pertain to the ring ge- 
ometry where r stands for the chord distance, r — 
{Na/ir) |sin(7ra;/A^a)|, x being the coordinate along the 
circumference. At large the electrons assume a 
Wigner-molecule configuration where they reside at the 
corners of a regular polygon. The effective low-energy 
Hamiltonian of such a state is given by 

^2 

H = —L"^ + j'^SjS.j + i +'^naflUJa, (3) 

j a 

where L is the centcr-of-mass angular momentum, Sj 
are electron spins, and are the occupation numbers of 
"molecular vibrations." (The same Hamiltonian applies 
to small- A'^ quantum dotsii^) At large Vg, the total mo- 
ment of inertia / and the vibrational frequencies are 
easy to compute because they are close to their classical 
values. Our task is to calculate J, which is more difficult. 

We show that the asymptotically exact relation exists 
between J and the PCF gQ{x) of a spin polarized ID 
system. For an ultrathin wire, C = \n{aB/R) ^ 1, it is 
particularly simple: 

^-^.^'(0), r.»i. (4) 

By virtue of Eq. the calculation of J reduces to an 
easier task of computing go{x). For large our final 
result has the form 

•^=7 :574 7 exp (-ryVZfl) , > 1, (5) 

where the values of r\ and k are given in Table n to- 
gether with the prediction of the FLA,-'^^-^'' for com- 
parison. [Klironomos et al^ independently obtained 
77 — 2.79805(5) for the case of a wire.] 

In the remainder of the paper we give the derivation 
of the above results and discuss their consequences for 
various experimental and theoretical questions. 

III. THREE-ELECTRON QUANTUM RING 

We start with the simplest nontrivial example: three 
electrons on a ring. Let < < 3a, j = 0,1,2 be their 
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coordinates along the circumference. We will compute 
the exchange coupling J between the j — and the j = 1 
electrons. In the FLA the third electron remains at rest 
while the other two interchange, as in a classical picture. 
However, we will show that the consistent calculation 
of J requires taking quantum fluctuations of this third 
electron into account. 

The system has the total of three degrees of freedom: 
the center of mass Xcm, the distance between the exchang- 
ing electrons x = xi — xq, and the distance between the 
third electron and the center of mass, X2 = X2 — Xcm ~ i- 
We can restrict the variables x and X2 to the fundamen- 
tal domain, \x\ < 3a/2, IX2I < a/2, where only the cho- 
sen two electrons can closely approach each other. This 
allows us to ignore exchanges involving the j ~ 2 elec- 
tron, and hence, its spin. Ignoring also the irrelevant 
center-of-mass motion, we obtain the Hamiltonian 



2/i 



-di- 



2M 



dl.^+Utotix,X2), 



(6) 



(7) 



where fi — 111/2 and M = 3/i. The potential term 

Utot = U{x) + Uo^t, 
C/cxt = C/[(3/2)(X2 + a)-x/2] 
+ U[{3/2){X2 + a) + x/2], 



has two minima in the fundamental domain, at x — ±a, 
X2 ~ 0. The minima are separated by a high poten- 
tial barrier at x = 0. They give rise to the two lowest- 
energy states of the system: the spin-singlet ground state, 
So+Si = 0, with an orbital wavefunction $s = ^s{x, X2) 
and a triplet with a wavefunction ■ Their energy split- 
ting is the desired exchange coupling J. In close analogy 
with Herring's classic treatment of H2-molecule^ (see 
also Ref. .22J for our = 3 "molecule" we obtain 



2fi,2 

J= — I dX2 



(8) 



where the (normalized to unity) "single-well" wavefunc- 
tion <I>i(x, ^2) is the ground-state of H3 with a modi- 
fied potential Utot Ui = [/tot(max{a;, 0}, ^2). Equa- 
tion is valid^^ to order O(J^); with the same accuracy. 



= -^[$i(x,A2)±$i(-a;,X2)]. (9) 

Let us discuss the form of the single- well wavefunction 
^i{x, X2). Near its maximum at a; = a, X2 = 0, it is a 
simple Gaussian, 



$1 {x, X2) (X exp 



where I = [hy^iUUa)] 



{x - af M 

~^ 2h 



21^ 



(10) 



1/4 



-1/4 



a is the amplitude of 



the zero-point motion in x (Utot written with a single 



argument is meant to be evaluated at X2 = 0). The 
quantity 



Lu{x) 



(^)] 



1/2 



(11) 



is real and positive in the classically forbidden region < 
X < a — I. Therefore, the tunneling barrier is the lowest 
and $1 is the largest along the line X2 = 0. Furthermore, 
it is easy to see that $1(2:, X2) rapidly decays at \X2\ > 
Since I <^ a, the following Gaussian approximation is 
justified in the entire fundamental domain of x: 



= (x) exp 



f,nix)Xl 



(12) 



Acting on this ansatz with Hj, and neglecting the sub- 
leading terms 0{l/a), we obtain the following equations 
on Q{x) and 4>{x): 



9?{x)-uj'^{x) 



\^^^dl-Utot{x) 



1/2 ' 

Kl{x) + E ) (t){x) 



[{2/^Ji)^Utot{x)] 
1 
2 



0, 



(13a) 
(13b) 



where /S.Utot{x) = Utot{x) — C/tot(a)- We wish to com- 
ment that the basic idea of the Gaussian ansatz has been 
used previously for computing exchange constants in '^He 
crystals)2& but the simple closed form of Eqs. H13|) has not 
been demonstrated. 

Comparing Eqs. ((Tn|l and (fT^ we see that near x = a 
function Mx) takes the form, 



4>{x) 



MQ (a) 



1/2 



(14) 



while fl{x) satisfies the boundary condition Q{a) = Lu{a). 
Resolving the 0/0 ambiguity in Eq. (|13a|l by L'Hopital's 
rule we furthermore obtain 



n' (a) = 



Lo'ia) 



1 



{(ill) i^ + u" («)/C/<^'.t(«)]} 



1/2- 



(15) 



Finally, using this formula and straightforward alge- 
braic manipulations one can fix the constant term E in 
Eq. ifTsKI) to be 



E = Utotia) + - [ujia) + luq] 



UJQ 



h 



— • (16) 



Next consider the interior of the classically forbidden re- 
gion, X <C a. Since the tunneling barrier A[/tot(a^) is 
large here, is a slow function of x. Taking advantage of 
the following expression for the PCF of a spin-polarized 
molecule, 



» N-l 

{x) = 2 n dA,<I>2(a;,A2,...,Ajv-i), |x| < 



3a 
Y' 

(17) 
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similar to Ref . (see also Appendix . Using such a 
representation, it is easy to show that 



(/)(O)/0'(O) ~ asIC, rs » 1, 



(24) 



which, combined with Eq. ((TH|l . yields Eq. |@J. With a 
bit more algebra, one can match Eq. I|22|l with another 
formula in region I, leading finally to Eq. (|3J) with r] and 
K given by 



FIG. 2: The PCF of a spin-polarized system (schematically). 
Regions I, II, and III are described in the main text. 



we find that Eq. JHJ entails the relation 



J 



0(0) 
4/i <t>'{Q) 



(18) 



Anticipating the discussion in Sec. II VI Eq. (|I7|) is written 
for an arbitrary N > 2, with the notation Xj = Xj — 
Xcm +(^^1^2j;')(a/2) being used; the PCF is normahzed 
as appropriate in the WC limit. 



3a/2 



go{x)dx = 1. 



(19) 



The dependence of go on x that results from this for- 
malism is sketched in Fig. |21 Near the x ^ a maximum 
(region III), go{x) is given by [cf. Eq. (^UJ] 



go{x) = (l/V^l) exp[-{x~af/f 



(20) 



In the region II the quasiclassical approximation applies. 
In particular, at <C x <C a the result for go{x) can be 
written in terms of the tunneling action, 

a 

Ssix) = \j dy [2^i^UUy)f^ , (21) 



and the appropriate prefactor, as follows: 



1 Q{a) hojQ 



27r n{x) U{x) 



1/2 



(22) 
(23) 



In comparison, the FLAii amounts to replacing fl{x) 
by a constant in the above equations, thereby effec- 
tively ignoring the quantum fluctuations of the dynam- 
ical variable X2. Finally, in an ultrathin wire, £ 3> 1, 
there is also region I, x < as, where the quasiclassi- 
cal approximation breaks down. Fortunately, at such 
x, Eq. H13b|) can be simplified, as il{x) ~ f2(0) and 
{/tot(a;) — U{x) + 2[/(3a/2). This enables us to express 
(j){x) and go{x) in terms of the Whittaker functionSf^ 



25/4 



dx 
a 



m. 



In{a) 



ea 



ea 



AC/tot(x) 



1/2 



C/tot(«) 



3/4 



(25) 
(26) 



Thus, for the N — i case we were able to reduce the orig- 
inal complicated three-body eigenvalue problem to rou- 
tine operations of solving an ordinary differential equa- 
tion H13a() and taking two quadratures, Eqs. (|23() and 
(I25|l . The resultant and k are listed in Table In 
comparison, the FLAii underestimates k, by about 50%. 
It gets 77 correctly but only for iV = 3, see more below. 

One more important comment is in order. The an- 
tisymmetry of the total fcrmion wavefunction imposes 
certain selection rules'^*' for the allowed values of L [see 
Eq. at a given total spin S. Thus, the lowest-energy 
L eigenstates for the two possible spin states of the = 3 
system, S = 1/2 and S — 3/2, are, respectively, \L\ = 1 
and 0. Since J ^ fi^ /I at large r^, the ground-state of 
the system is the L — spin-quartet, in agreement with 
prior numerical work^^iS 



IV. CASE iV > 3 

In a system of more than three electrons, the Hamil- 
tonian that governs the important degrees of freedom, x 
and X = (A2, . . . ,Xm-i)\ becomes 



Hn — 
where 



2/i 2 



(M-V2ax)^(M-V2ax)+c/tot, 

(27) 



-1/2 



1 %-(l-v/27]v) 



N -2 



(28) 



Again, the potential energy has two minima separated 
by a large barrier. The single-well function <I>i(a;,X) can 
be sought in the form 



$1 — (j){x) exp 



-4rAXtMi/2f^(a;)Mi/2AX 
2h ^ ' 



(29) 



where AX = X — X*. In the language of the quan- 
tum tunneling theory, X*(a;) is the instanton trajectory 
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and fl{x) is a matrix that controls Gaussian fluctuations 
around the instanton. Our goal is to compute them. The 
usual route is to parametrize the dependence of X* on 
X in terms of an "imaginary-time" r, in which case x(t) 
and X* (r) must minimize the action functional 



could expand At/tot in Eq. (|30|l to the second order in 



^{drxr 



(30) 

subject to the boundary conditions a;(0) = 0, a;(oo) = a, 
and X(oo) — 0. Henceforth [/tot is always meant to be 
evaluated on the instanton trajectory and AC/tot stands 
for the difference of its values at a given t and at r = oo. 
Repeating the steps of the derivation for the = 3 
case, we find that Eq. (|13a|l for Jl is still valid once 
we define u; to be a positive-definite matrix such that 
(jj^ = M~^/^HM~^/^, where H is the matrix of the sec- 
ond derivatives S^j — dxidxjUtot- This equation has an 
equivalent but more elegant form in terms of r: 



drVt = 02(t) -u;2(t) 



(31) 



Also, for practical calculations, it is convenient to take 
advantage of a formula 



1 r)2 1 

(32) 

As for Eqs. I)13b|) and H23|l . they require only the re- 
placement 57 tr r2, 

|^52-C/tot(x)-^trf2(x)-f i?|0(x) =0, (33) 

while in Eq. H26|) one has to replace 17 with det f2. 

Finally, instead of Eq. (^5)) we have -q — 2SN/y/2r^, 
consistent with the notion that the main exponential de- 
pendence of J is always determined by the tunneling ac- 
tion. 



V. CALCULATION OF THE INSTANTON 

A few properties of the instanton follow from general 
considerations. The dimensional analysis of Eq. H30() 
yields Sn oc so that is indeed just a constant. 

Also, from the symmetry of the problem, Xjv+i-j(T) = 
—Xj{T). Thus, in the special case of = 3, the instan- 
ton trajectory is trivial: X2 = 0, i.e., the j = 2 electron 
does not move. This is why we were able to compute ^3 
in a closed form, Eq. (|21|) . For > 3 the situation is 
quite different: all electrons [except j — {N + l)/2 for 
odd N] do move. In order to investigate how important 
the motion of electrons distant from the j = 0, 1-pair is 
let us consider the N = 00 (quantum wire) case, where 
the far-field effects are the largest. If Xj^s were small, we 



Xj to obtain the harmonic action 



1 m f dk f duj . 



-^p(fc)], (34) 



where Ukui is the Fourier transform of Uj{T) = Xj — x^, 
electron displacement from the classical equilibrium po- 
sition = (j — l/2)a, i G Z, and 



4.15 



So 




(35) 



is the plasmon dispersion in the ID WC (the logarithmic 
term is due to the long-range nature of the Coulomb in- 
teraction) . Minimization of Sh with the specified bound- 
ary conditions yields 



Uj{t) cx 



(a;°)2 -I- 



^InMll 
2 



(36) 



Substituting this formula into Eq. H34() . we find that the 
contributions of distant electrons to Sh rapidly decay 
with \j\. Since Uj (r) 's are small at large j and t we expect 
that these j- and r-dependencies are rendered correctly 
by the harmonic approximation. Thus, a fast conver- 
gence of 77 to its thermodynamic limit is expected as N 
increases. 

Encouraged by this conclusion, we undertook a di- 
rect numerical minimization of S for the set of N listed 
m Table n using standard algorithms of a popular soft- 
ware package MATLAB. The optimal trajectories that 
we found for the case A^ = 8 are shown in Fig. O In 
agreement with our earlier statement, Uj{0) reach a fi- 
nite fraction of a. This collective electron motion lowers 
the effective tunneling barrier and causes 77 to drop be- 
low its FLA value, although only by 0.7%, see Table 
Further decrease of 77 as A" increases past 8 is beyond the 
accuracy of the employed minimization procedure. 

Let us now discuss the prefactor k. In the inset of 
Fig. U] we plot trr2(a;) computed by solving Eq. H31() nu- 
merically. To reduce the calculational burden, we set 
X*(t) — > instead of using the true instanton trajectory. 
The error in k incurred thereby is ~ 2% (see Sec. IVII|I . 
In comparison, the FLA, where trrj(x) = const, yields 
K about 50% smaller than the correct result, similar to 
A^ = 3. 



VI. ALTERNATIVE DERIVATION OF THE 
EXCHANGE CONSTANT 

The relation (0J between the exchange constant J 
and the PCF of the spinless system — which is one 
of our main results — was derived in Sec. IIIII assum- 
ing ^ 1. In this section we show that in ultrathin 
wires, C = ln(aB/i?) 3> 1, this relation holds under a 
more general condition: with an accuracy 0{1/C), it re- 
mains valid as long as the strong spin-charge separation 
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FIG. 3: (Color online) The instanton trajectories of 1 < j < 4 
electrons in the N = 8 Wigner molecule on a ring. Inset: 
trf2(a;). The units of r and CI are ^/2a/so and its inverse. 



persists, J ^ Ep ~ fi?/ma?. In particular, it applies 
at Tg ~ 1, i.e., in a parametric regime which has been 
notoriously difficult for a controlled theoretical analysis. 
Therefore, we think that Eq. Q amounts to some def- 
inite progress. The fundamental reason for the broad 
domain of validity of Eq. Q is that in ID the 1/r fall- 
off of the Coulomb potential jSJ constitutes a sufficiently 
rapid decay. In fact, the relations of type of Eq. are 
generic for ID models in which electrons interact via the 
potential U(r) ~ 1/r" with a > 1. Known to us physi- 
cal realizations of a > 1 cases include (i) a = 3, which 
occurs when the bare Coulomb interaction is screened 
by a nearby metallic plane, and (ii) a = 2, where the 
screening is accomplished by a half-plane, see Sec. IVIII 
below. 

In order to derive Eq. Q for the unscreened Coulomb 
interaction [a — 1) we make a connection with a previous 
work of one of us^^ where it was shown that in ultrathin 
wires an unusual correlated regime — Coulomb Tonks gas 
(CTG) — exists in the window l/£ ^ Ts ^ 1. In the 
CTG, unlike in a WC, the long-range tails of the Coulomb 
potential have no important effect on the spin degree of 
freedom. To the leading order in r^, the Coulomb inter- 
action Q acts simply as a strong short-range repulsion. 
This allows one to characterize its effect solely in terms 
of the transmission coefficient t[q) for a two-electron col- 
lision with the relative momentum q. Based on this ob- 
servation, the following formula for J was derivedSS 



f dq 

J =— —ga{q)qlnit{q), J Ef 
m J In 



(37) 



To see how it leads to Eq. we first note that for the 
potential |(2J) the transmission coefficient has the form— 



iq 



iq — {m/h?)c{q) ' 



9> 



1 

as 



(38) 



where c{q) = 2(e^/e) In (l/2gi?). The slow (in this case, 
logarithmic) dependence of c on g is crucial here and 
it is precisely due to the aforementioned fast decay of 
U (r) with distance. As explained in Ref. 22, in the CTG 



regime the integral in Eq. H37|) is dominated hy q ^ 1/a; 
therefore, to the order 0{\/C) one can replace c{q) by 
c(l/a) and then by 2(e^/e)£. Finally, it is permissible 
to neglect iq in the denominator of Eq. H38|) . thereupon 
Eq. l|H7|) becomes identical to Eq. igj. 

Our next task is to show that Eq. holds not only 
in the two limiting case (CTG and WC) but also every- 
where in between. We do so with the help of a formal- 
ism that can be considered a continuum analog of the 
Schrieffer- Wolfe transformation in the Hubbard model, 
in particular, in its ID version;-^; 

Let J G Z label electrons in the order of increasing in- 
stantaneous coordinates Xj. Then, by construction, the 
coupling constant J refers to the exchanges of nearest- 
neighbors, say, i = Q and j = 1. This coupling is sup- 
pressed because in order to exchange, the two electrons 
have to approach each other very closely and at that 
moment their mutual Coulomb repulsion U{x) is very 
strong. Let ac be a suitable short-range cutoff. Since 
we are interested in the physics at the energy scale J, it 
should be possible in principle to "integrate out" the pro- 
cesses at the much higher energy scale U{ac)- In doing so 
one would trade the bare Hamiltonian H for a renormal- 
ized one, H + , which however yields the same value of 
the effective coupling J. Traditionally, this renormaliza- 
tion program is accomplished by incrementally increasing 
Uc from zero to a, the natural upper limit. In particular, 
the Hamiltonian can be viewed as the final result of 
such a renormalization. However, it is fully legitimate to 
stop the process while Oc is still much smaller than a. At 
this stage the charge correlations at scales shorter than 
ttc can not longer be studied because the corresponding 
degrees of freedom are already integrated out. This im- 
plies that the form of is not unique: iJo- may differ in 
their structure on scales x < ac as long as the resultant J 
is unchanged. Therefore, it should be possible to choose 
Hrr in the form 



= [/„(a;)(SoSi - 1/4) + AcSix). 



(39) 



It is especially convenient to take the limit Ac — > oo, 
in which electrons behave as impenetrable particles. Al- 
though thus renormalized Hamiltonian forbids the pair- 
exchange in the real space, if the choose Ua {x) correctly, 
the spin exchange rate J, which is the only observable 
manifestation of an interchange of identical particles, 
would be preserved. A familiar example of such a 
on a lattice is the added term 



1 



(40) 



complemented with the no-double-occupancy constraint 
(i.e., U ^ U + AU = 00) in the theory of a large-C/ ID 
Hubbard model. 

The choice of Ucr{x) in Eq. is again not unique. To 
readily get the result we want, we require that (i) Ua{x) 
decays exponentially with a; and (ii) its effect on any 
wavefunction $(a;) that vanishes at a; = is adequately 
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described by the first-order perturbation theory. Under 
such conditions the spin dynamics is captured correctly 
to the leading order in 1 /£ as long as Ua {x) satisfies the 
constraint 



dx x'^Ucr{x) 



(41) 



The proof is relegated to Appendix ^ From this point 
the derivation of Eq. is straightforward. After the 
renormalization, the infinitely strong ^-function included 
in causes the orbital wavefunction of the system to 
vanish at x = 0, so that the perturbation theory in Ua 
now applies. To the leading order we can compute the 
sought exchange constant J by averaging Ua over such a 
wavefunction or equivalently over the renormalized PCF 

9rcn {x^ . 



which is similar to the celebrated Calogero-Sutherland- 
Moser (CSM) modelf^ except for the short-range cutoff 
R in U{r). Without this cutoff (i.e., for i? = 0) the 
potential U is impenetrable, so that J — 0. With a fi- 
nite cutoff a weak tunneling through the potential barrier 
U (r) is possible, which imparts the system with a spin 
dynamics and has interesting experimental implications 
(see Appendix lBjl . But to finish with the theoretical part, 
let us first focus on the original CSM model, R — 0. This 
is a good test case because, on the one hand, a number of 
exact analytical results are available here. On the other 
hand, the model also possesses a WC regime, at A » 1, 
so that the corresponding PCF go{x) should be calcula- 
ble by our method. According to the exact resultS j^^i?^ 
the PCF behaves as 



9oix) 



2A 



0. 



(44) 



J = -j^ I Ua{x)g,cn{x)dx. 



(42) 



(The factor 1/2 accounts for the two-body nature of the 
interaction.) Let us discuss grenix) in more detail. First 
of all, under the conditions of a strong spin-charge sep- 
aration the PCF a spinful (g) and spinless (50) systems 
are always very similar. Before the renormalization (for 
the original Hamiltonian) both functions are peaked at 
x ~ a and rapidly decrease as a; ^ 0, albeit only go{x) 
vanishes at x = 0, while g{Q) is simply very small. Once 
Ha is added, g{x) gets replaced by the function gi-cn{x), 
which vanishes at x = and thus is nearly identical to 
the original go x. The PCF goix) of the spin- 

less system does not change appreciably. Since Ua{x) 
is short-range, the integral in Eq. (|42|l is dominated by 
small X. Therefore, we can substitute the expansion 
grcn{x) ~ goix) ~ (^g (0)x^/2 in place of the full gron(a^)- 
Comparing the result with Eq. (|41|l . we arrive at Eq. 

This alternative route to the derivation of Eq. has 
the advantage of being valid in a broader range of Vs 
compared to Herring's method2& originally designed for 
the quasiclassical case Vs ^ I. It may also be more 
intuitive because the Hubbard-model expression H4U|I is 
widely known. Finally, it can be easily generalized to 
other interaction types, as discussed in the next section. 



VII. AN INTEGRABLE MODEL WITH A 
SCREENED COULOMB INTERACTION 

Our instanton calculation in Sec. II VI is expected to be 
asymptotically exact in the limit of large Vs . However it 
is a good practice to verify whenever possible that any 
given calculation is free of inadvertent arithmetic or cod- 
ing mistakes. With this in mind we applied our method 
also to the interaction law 



Uir) 



A(A - 1) 

m (r + R)^ ' 



(43) 



For the three-electron molecule on a ring the coefficients 
C and m: can be computed in a straightforward manner 
directly from the exac^ three-body wavefunction: 



A^ = 3, A> 1 : C = 



8tt 

9^3' 



(45) 



For N — 00 the derivation is much more involved, but 
the result^iSi is also known: 



C 
2tt\ 



r3(A + i)v/3M 



-I 1/(2A) 



r(2A + i)r(3A + 1) 



where T{z) is Euler gamma- funct ion In the WC limit 
the first formula reduces to 



C 



3V3' 



A> 1. 



(47) 



The calculation of go{x) by the instanton method is virtu- 
ally the same as in Sec. Ivl except for two minor changes. 
First, the quasiclassical approximation is valid down to 
X = 0, i.e., there is no region I in Fig. [21 Second, to han- 
dle the logarithmic divergence of the tunneling action, 
inherent to the CSM model, one has to impose a differ- 
ent boundary condition for the instanton, x(t = 0) = z, 
which makes Sn a function of z [similar to S3 being a 
function of x in Eq. 121|) ]. At small z it behaves as 



Sn ^ y/X{X - 1) [ln(a/z) - AS"]. 



(48) 



Replacing ^A(A - 1) by A - 1/2, we find [cf. Eqs. ^ 
and ^] 



C 



MO)-AS 



iVX/c 



I det fl{a) 
detf2(0) 



(49) 



(it is easy to see that x is a, A-independent constant) . 
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Remarkably, for N = 3 the solution of Eq. (|13a|l and 
all the necessary quadratures can be done analytically. 
Using the expression for total potential energy 



f/tot (x) 



A(A - 1) 

9 m 

1 



sin (7ra;/3a) cos^ (7rx/6a) 



(50) 



we derive (for A » 1) 



0(x) 



■ . /7ra:\ / nx 



3a) 



(51) 



3 ma^ cos^ {Trx/6a) ' 



which allows one to recover the exact result Eq. H45|l . 

To check if we can also reproduce Eq. (|47(l for N = oo, 
we calculated AS" and x numerically for 4 < iV < 10 and 
fitted them to cubic polynomials in 1/iV. Extrapolating 
the fits to iV = oo, we obtained exp(AS') = 1.6438(4) 
and X = 1.04 compared to the exact values 1.6434 and 
1.02, respectively. We attribute the 2% discrepancy in k 
to our choice not to use the true instanton trajectory in 
Eq. H31|l . Apparently, our method has successfully passed 
this test; thus, our results for the unscreened Coulomb 
interaction (Table H)) should also be reliable. 

Let us now consider the case with a small but non- 
zero cutoff, < i? <C a. As explained in Appendix 
this model can be relevant for the existing experimental 
setup of Auslaender et al^ provided the electron density 
can be made low enough. Therefore, the calculation of 
J has an independent significance. This calculation is 
again similar to that for the case of the bare Coulomb 
interaction. Sec. I VII except that instead of Eq. H41|) the 
constraint on the effective potential becomes 

J2A ^ / dxx^'U, (x) - -^^^i?2A-i. (52) 
Jo A - 1 TO 

(For derivation, see Appendix 1X1 ) Consequently, Eq. 
is replaced with 

J = hxa^'^^goiac), i? < < a. (53) 

Since R is parametrically smaller than a, it is permissible 
to substitute Eqs. 144(1 and H46|l into the last formula, 
which yields 



J 



1 



R 



2A-1 



1 (2^A)'^+^r3(A+l) 
2A(A-1) r(2A + l)r(3A+ 1) ■ 



(54) 
(55) 



From Eq. l(51|l one concludes that the strong spin-charge 
separation in the model in hand may arise for an arbi- 
trary A > 1 provided R/a <C min{l, A — 1}. In the ex- 
perimental setup of Auslaender et a/. where this model 
can be realized physically, we estimate R ~ O.Ta^ and 



A = 1.5-2, see Appendix FEl so that 7 = 36-40. Our for- 
mula predicts that as the electron density is reduced, a 
fairly rapid algebraic fall-off Eq. H54|l of J should become 
observable, provided the disorder effects do not intervene. 
The crossover from high to low-density behavior is dis- 
cussed in more detail in the next section where we give 
arguments that the strong spin-charge separation sets in 
already at rather modest r^. 



VIII. VALIDITY DOMAIN OF THE OBTAINED 
FORMULAS AND NUMERICAL ESTIMATES 

A. Intermediate densities 

Returning to the case of the unscreened Coulomb inter- 
action, an important practical question is the relevance 
of the obtained formulas for a more common experimen- 
tal situation of Ts ~ 1. Although Eq. Q reduces this 
question to an easier problem of calculating the PCF of 
spinless fermions, at such it still has to be done numer- 
ically. Deferring this effort for a possible future work, we 
attempted to see whether one can estimate J at Ts ~ 1 
by a simple interpolation between the high-density limit 
Ts ^ 1, corresponding to the CTG,^^ and the low-density 
WC limit, Ts 3> 1, studied in this work. 

Using Eq. 13711 w ith go{q) computed to the linear order 
in Ts (per Ref. l23) we obtained for the CTG 



eas 24rf£ 



2(3W) 
3^2 



O 



(56) 



which is valid for l/£ ^ j's ^ 1- We plotted this depen- 
dence in Fig. Altogether with the prediction of Eq. (jSJ for 
the infinite-wire WC. As one can see, a smooth interpo- 
lation between the two is entirely possible. In fact, the 
curves match almost seamlessly around ~ 0.9 where a 
tiny gap has been left on purpose to separate them from 
each other. However, we must warn the reader that at 
such Vs the high-density approximation of Eq. H56I) cer- 
tainly needs higher-order corrections (taken literally, it 
would give negative J at > 1.15). Similarly, Eq. is 
likely to have its own non-negligible corrections at such 
Tg. Thus, one should regard Fig. 0] as the indication that 
the crossover from the CTG to the WC occurs at — 1- 
2. To put it another way, we expect that our Eq. ^ be- 
comes quantitatively accurate at Ts > 2. This conclusion 
seems to be in agreement with the numerical result of 
Ref. that the ID WC is well-formed already at > 4. 
To get an idea of the size of J, we can use parameters 
Ts =4, as — 1.5 nm, and e = 1, which are relevant for 
the experiment of Jarillo-Herrero et al^ For this set of 
numbers Eq. (0 yields 



J w 1 K. 



(57) 



Unfortunately, the lowest temperature in that experi- 
ment was 0.3 K, so the exchange correlations may have 
been strongly affected. We leave the investigation of the 
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/ 0.2 




FIG. 4: Exchange constant in units of /eas as a function 
of Ts for as/R = 100. The left curve is computed according 
to Eq. H56|l ■ the right one is per Eq. lO. 



finite-temperature effects for a future study. We also 
hope that lower temperatures can be achieved in the next 
round of experiments, so that we would be in a better po- 
sition to check our predictions. 



B. Thicker wires 

Another question that we wish to briefly address is the 
behavior of the exchange constant in a wire of radius 
o,B- This parametric regime is more common than 
that of ultrathin wires, R <^ as, and so it needs to be 
investigated. 

In regards to the WC limit, the correction to param- 
eter r] of Eq. (O due to a finite R has been studied by 
Klironomos et al^ These authors pointed out that in 
the wires with radius R ^ as the problem of spin ex- 
change cannot be treated as a purely ID one from the 
outset. They found that the tunneling trajectories of 
the primary electron pair avoid the head-on collision by 
deviating off the wire center in the transverse direction 
by a certain amount rg each. For the parabolic confine- 
ment potential Ucon{r±) = {h'^/2mR^){r±/R)'^ it is easy 
to find the estimate Tq ~ {R'^/asY^^ by balancing the 
Coulomb repulsion and the lateral confinement forces at 
the point of the closest approach. As long as rg ^ a, it 
remains possible to modify the interaction potential U{x) 
to take this effect into account. [At the simplest level, it 
is sufficient to replace R by 2ro in Eq. (0).] Thereafter, 
one can still treat the problem as IDi^ 

In thick wires the relationship between the PCF and 
the exchange constant would differ from Eq. I^J. Of pri- 
mary importance here is the small-distance behavior of 
U{x). When ro S> a^, the interaction potential becomes 
so smooth near x — Q that it can be approximated by a 
constant J7(0) ^ e^/2ero. Following the method outlined 
in Appendix fXl one arrives at the relation 



We see that the factor C in Eq. Q gets replaced with 
as/ao « lia 

In the case of intermediate-width wires, R ^ as, one 
can use uq = as for order-of-magnitude estimates; how- 
ever, the exact value of the coefficient oq and therefore k 
in Eq. lO would depend on the details of the confinement 
potential ?7con(r±)- 



IX. CONCLUSIONS 

In this paper we have developed a method for a con- 
trolled calculation of the exchange constant J in corre- 
lated ID electron systems to the leading order in the 
large parameter r^. This method is able to reproduce 
a number of exact results for integrable models, which 
is a strong indicator of its validity. In comparison to a 
prior uncontrolled calculation^^ of the same quantity for 
the Coulomb case, our method brings a small correction 
of about 0.7% to the coefficient in the exponential term; 
however, the correction to the pre-exponential factor is 
large, about 50%. 

Our another main result — the relation between J and 
the pair-correlation function of spinless electrons — is a 
promising ground for an attack on the regime of interme- 
diate concentrations, Ts ~ 1, which has so far been diffi- 
cult for a controlled analytical study. We demonstrated 
that a very smooth interpolation between the high and 
low-density asymptotic is possible based already on the 
results in hand (Fig. . 

Finally, we have applied our method to a number of 
realistic interaction potentials corresponding to different 
geometries of experimental setup. Apart from nearly 
impenetrable Coulomb interaction in ultrathin quantum 
wires and quantum rings, which has been the primary 
focus of this work, we have studied penetrable smooth 
interactions that take place in the wires of intermediate 
width and a screened interaction that is realized, 

e.g., when the wire is located near a metallic half-plane. 

These findings have direct implications for the en- 
ergy spectroscopy of spin-splitting in nanoscale quan- 
tum rings^i^ and ID quantum dots.^'"'^°'^^ Some esti- 
mates were given in Sees. IVIII and IVIIII above. They 
should be especially accurate in the ultrathin- wire limit — 
C — ln{aB/R) S> 1 that can be achieved in carbon- 
nanotube field-effect transistors with high-A: dielectricsi^ 

In long ID wires, the exchange coupling J determines 
the velocity of the elementary spin excitations 



TT Ja 
2T' 



(59) 



J 



2e 



aog'oix), ao = 



VW7(0) 



(R^aB)'/'. (58) 



which can nowadays be measured by tunneling, angle- 
resolved photoemissionji^ or deduced from the enhance- 
ment of the spin susceptibility and electron specific 
heat.^^ Our result for reads (cf. Tabled 

w^/i^F = IV.Sr^/^e^"^//:, 77 = 2.7978(2), (60) 
where vp — {t: /2){h/ma) is the Fermi velocity. 
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Wavefunctions $s and $t are the linear combinations 
of u and v that are even and odd in x, respectively. For 
X at which Eq. ljAl|l is valid, — — const x v{x). 
Using this property and some trivial algebra, we get 



dxUa{x)v'^{x) 



(A4) 



APPENDIX A: SPIN-DEPENDENT PART OF 
THE EFFECTIVE POTENTIAL 

In this Appendix we derive Eqs. l(IT|) . and (|S5|l . 

Let X = xi — xo he the relative distance between the two 
electrons for which we want to compute the exchange cou- 
pling. Consider the orbital part of the many-body wave- 
function <I>, which is a function of x and A^— 1 other coordi- 
nates. In principle, $ is different for different many-body 
spin states. Nonetheless, under the condition J ^ Ep 
where only nearest-neighbor exchanges are important, at 
small X function <I> depends primarily on the total spin 
of the given pair. The triplet-state wavefunction $t(x) 
is odd and therefore $t(0) = 0. In the singlet state we 
have ^s{x) which is even and finite at x = (cf. Sec. lIIip . 
However, because of the strong Coulomb repulsion <i>s(0) 
is very small. At a; <C a the two electrons of interest are 
under a large Coulomb barrier and it is legitimate to as- 
sume that $t and $s are functions of only x, the fastest 
variable in the problem. Our goal is to choose 11^ (x) that 
reproduces the difference between $s and $t at the level 
of the first-order perturbation theory. This is achieved if 
there exists a window Oc ^ a; ^ a, where the following 
equation is satisfied: 

oo 

- $,(a;) = J dx'Gix,x')U„ix')^tix'). (Al) 



Here the Green's function G is the solution of 

fi2 



G{x,x')^-6{x-x'), (A2) 



where E ~ e^/ea is the total energy. At the end of the 
calculation one can verify that the final result [Eq. (|41|l ] 
holds independently of the precise value of E up to terms 

0(1//:). 

Our next step is to use the usual representation of G, 



Gix,x') ^ - 



2/i J u{x)v{x'), x > x' , 
v{x)u{x'), x < x' , 



(A3) 



For the Coulomb case [Eq. (O] u and v can be expressed 
in terms of Whittaker functions^^ W^fz) and Ma{z) (see 
an earlier remark in Sec. IIIII and Ref. l23): 



(A5) 



u{x) = PF_j,.,i/2[«(a; - R)/b], 
v{x) = M_i^^i/2[i{x - R)/b] 

-M-«.4/2(-*fi/fe)"(-^)/"(0), 



where b = ?i/(8/xi5)^/^ and = b/as ^ ^rj > 1. 
Using the relations^S l^-ii/,i/2(0) — l/r(l + ii') and 
M_i^ i/2{iz) ~ at z <C 1, we find Q = — w)b. 

Substituting these results into Eq. (|A4I) and dropping 
terms that are small in the parameter R/as <C 1, we 
recover Eq. l|lT|l . 

Similarly, for the CSM model [Eq. u and v are 

given by 



u{x) = ^/xH[^}^^2{qx), 
v{x) = VxJx-i/2{qx), 



(A6) 



where q = y/mE/h, and Hi, , are the Hankel and the 
Bessel functions of the first kind,^^ respectively. Their 
Wronskian is given by u} = l/i7r(A — 1/2), and 

their known asymptotic behavior leads in this case to 
the formulas 



u{R) ~ r{X-l/2)R^-^/iTT, 
u\R) ~ (1 - \)u{R)lR, 
w(i?)~i?Vr(A + l/2), 



(A7) 



so that the constraint on Ucr{x) is expressed by Eq. H52|l . 

Finally, in thick ID wires, where U{x) ~ const 
at small x (Sec. IVIIIjl . the solutions of Eq. ljA2|l are 
u{x) = exp(— x/flo) and v{x) = sinh(a;/ao), where uq = 
h/ y^mU{0). These formulas leads to Eq. (jJHJ • 



APPENDIX B: PHYSICAL REALIZATION OF 
THE CALOGERO-SUTHERLAND-MOSER 
MODEL 



in terms of two linearly-independent solutions of the (ho- 
mogeneous) Schrodinger equation (|A2p . The solution 
u{x) initially decays exponentially with x and then turns 
into a wave propagating towards x = +oo. The other so- 
lution, v{x) has a node at a; = 0, exhibits an exponential 
rise and finally becomes a standing wave. The quantity 
Q = uv' — u'v is their Wronskian. 



The CSM interaction law can be realized experimen- 
tally if the ID wire is positioned nearby a metallic half- 
plane, parallel to its edge, see Fig. [3 In this configu- 
ration the interaction potential has the Coulomb form 
[Eq. ||2Jl] only at small x. Below we show that at larger 
distances, |x| ^ D, the interaction changes to the CSM 
law U{x) ~ A/x'^, due to the screening effect of the 
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metal. For simplicity, let us assume that the metallic 
half-plane is an ideal conductor and that x ^ R. In this 
case U{x) = V{x,D,9), where V is the solution of the 
following electrostatic problem: 



VV(r) 



1 92 



\dx'^ 



p dp^ dp 



eD 

V{r) =0, = ±7r 



p2 dip^ 

5{x)5{p ~ D)5{ip ~ 9 



V{x,p, (fi) 



(Bl) 
(B2) 



Using standard methods (cf. Ref. 36 and references 
therein) one can express V in terms of elliptic integrals. 
However, for a general 9, such expressions are not very 
illuminating. An exception is the case 9^0, where V 
takes the form 



V{r) = 



ne r 



D 



^2Dp{l + cosip) 
arctan-i^ j — , (B3) 



where D = (0, 0,D). This formula can be verified by 
the direct substitution into Eqs. (jBljl and ljB2|l . For the 
intra-wire interaction law we obtain 



N 2 e2 2D 
U [X) — — - — r arctan - — - , 
7re \x\ \x\ 



(B4) 



so that A = 4De^ /ire in this case. We can also show that 
for an arbitrary — tt < 9 < tt the coefficient A is given by 



First we note that the Fourier transform of U must have 
the form U{q) ~ U{0) — TTA\q\ at small q. On the other 
hand, U is given by the generalized Fourier series 



9^2 °° 
n— 1 



i-ircosn9]Ir,/2{qD)K,,f2{qD), 



where Ii, and are the modified Bessel functions of 
the first and the second kind, respectively. The non- 
analytic |g|-contribution comes only from the n = 1 
term in the series. Using the formulasSS Ii/2iz) — 
(2/7rz)^/2 sinhz and isri/2(z) = (7r/2z)^/2 exp(—z), we ar- 
rive at Eq. (|B5|I . 

The geometry of Fig. Elb) has been realized in the 
experiment of Auslaender et al.^, where the role of the 
metallic half-plane was played by a two-dimensional elec- 
tron gas. In contrast to our idealized electrostatic model, 
the two-dimensional metal has a non-zero Thomas-Fermi 
screening radius ttf- We expect that in this situation 
the only change in Eq. ljB5|) is the replacement of D by 
D + rxF • This yields the following estimate for the pa- 
rameter A in Eq. H43|l : 



A = 



1 2 D -f ttf 

7 + 

4 TT as 



(1 + cost 



(B6) 



A 



2D 



(1-f cos^). 



(B5) 



In the experimental setup of Ref. \^ D ^ as ^ ttf , so 
that A = 1.5-2. 
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Significant changes in the calculation would be required 
only in very thick wires with radius R > 11 asi^ However, 
in such wires the ID WG limit can be achieved only at > 
20 where the exchange coupling J would be unmeasurably 
small, see Eq. 0. At smaller r^, other physics is likely to 
come into play, e.g., the ID WG may undergo a transition 
into a multi-row WO—*— or into a liquid with multiple 
subband occupation. 

It is interesting that in the opposite limit of a short-range 
U{x), e.g., U(x) — cS{x), the relation between J and go 
is formally similar, with ao = 2h? /mc, see Eqs. 13711 and 



